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For heavy-ion collisions at RHIC a scaling behavior is found in the dependencies on azimuthal 
angle <j> and impact parameter b for pion production at high pr essentially independent of the 
hadronization process. The scaling variable is in terms of a dynamical path length £ that takes 
into account detailed properties of geometry, medium density and probability of hard scattering. It 
is shown in the recombination model how the nuclear modification factor depends on the average 
£(4>, b). The data for ir° production at px = 4-5 and 7-8 GeV/c at RHIC are shown to exhibit the 
same scaling behavior as found in the model calculation. Extension to back-to-back dijet production 
has been carried out, showing the existence of £ scaling also in the away-side yield per trigger. At 
LHC the hard-parton density can be high enough to realize the likelihood of recombination of shower 
partons arising from neighboring jets. It is shown that such 2-jet recombination can cause strong 
violation of f scaling. Furthermore, the large value of Raa that exceeds 1 can become a striking 
signature of such a hadronization process at high energy. 

PACS numbers: 25.75.-q, 25.75.Gz, 24.85. +p 
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I. INTRODUCTION 

Recent analysis of the high-statistics data on ir° pro- 
duction in heavy-ion collisions has provided a wealth of 
information on the nuclear modification factor (Raa) as a 
function of centrality (c), azimuthal angle (<p), and trans- 
verse momentum (pt) The complexity revealed is 
clearly in need of some organizational simplification, if 
useful insight is to be gained from it to advance a theoret- 
ical understanding of the various features. For example, 
the pt measured ranges from 1 to 10 GeV/c; the con- 
ventional wisdom is to separate the high and low regions 
and treat them with different dynamics. The <fi depen- 
dence is studied by determining the second harmonics 
(1)2), whose dependencies on c and pr are presented as if 
they are variables orthogonal to <ft, although it is known 
that c and <f> are correlated by geometry and px depends 
on path length. On the theory side treatments of funda- 
mental issues at the microscopic level cannot incorporate 
experimental details without sacrificing clarity and pre- 
cision. In the interest of finding simplifying features of 
the data we examine in detail the geometrical aspect of 
the problem associated with the propagation of a hard 
parton through the medium and explore the possibility 
of scaling behavior in a variable, £, that involves both im- 
pact parameter (b) and <f>. On the basis of that result we 
then organize the dynamical problem of hadronization in 
a way that can exhibit the scaling behavior. With that 
behavior derived on a theoretical basis in certain approxi- 
mation, the data on Raa can then be shown in the same 
format to check the relevance of the model-dependent 
finding. A form of universality does exist in reality. 

At low pt where thermal and semihard partons dom- 
inate, a study of the surface factor S(c/>, b) that relates b 
and (j) has led to the elucidation of the <f> dependencies 
of RAA(4>,b) and ridge yield Y R (4>,b) Q. At higher pr 



shower partons become more important Q . Accordingly, 
the path length of a hard parton in the dense medium 
becomes the major factor that affects Raa- That length 
depends on <j>, b, and the point of creation of the hard 
parton, the probability of which depends on the nuclear 
overlap function. All those quantities are calculable from 
geometrical considerations. Of course, at some point dy- 
namics enters the picture. Recombination can handle 
well the late-stage hadronization process at all pr 0]. 
Our primary objective is to expose the simplifying fea- 
tures of the nuclear complexity that a hard parton expe- 
riences at the early stage of the production process. It is 
complementary to the study of the jet quenching process 
in QCD that assumes a simple path length L [4|, |5| . The 
two meet at the point where momentum degradation is 
expressed in terms of the scaling variable £((/>, b). The 
connection is complicated by the role that thermal par- 
tons play in the hadronization of shower partons at the 
intermediate px range. Nevertheless, a simple behavior 
of Raa in terms of £ can be determined theoretically and 
experimentally. 

If indeed a universal behavior of Raa can be found 
in its dependence on a dynamical path length, then it 
should be natural to follow up the investigation with a 
study of the properties of back-to-back dijets, since the 
dependence of the yield on the momenta of the trigger 
and away-side particles (pt and p&, respectively) is closely 
related to the path lengths traveled by the hard parton 
and its recoil. A detailed study of the dependence on 
p t and pi, has been carried out earlier, but with <f> av- 
eraged over Q. Now, with specific 4> values taken into 
account, the joint dependence on b and <fi presents com- 
plications that must be simplified, if useful comparison 
between theory and experiment is to be facilitated. 

Dijets mentioned above are directed in opposite direc- 
tions, but two jets going in the same direction would be 
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of tremendous importance when collision energy is very 
high as at LHC. Pion production arising from 2-jet re- 
combination should overwhelm 1-jet fragmentation when 
the rate of hard scattering is significantly increased at 
high initial energy. Our procedure of studying the <f> and 
b dependence can readily be extended to the coalescence 
of shower partons from parallel jets. Our result shows 
the possibility of dramatic departure from the usual ex- 
pectation on the nuclear modification factor. 

In this work we shall focus on the production of pi- 
ons, since the high-statistics data pj that can check our 
analysis are for ir° only. The extension to the production 
of protons and other hadrons can be done in the same 
framework, but will not be carried out here. 



II. INCLUSIVE PION DISTRIBUTION AT 
FIXED AZIMUTHAL ANGLE 



We shall work in a pr range where thermal-thermal re- 
combination is a small correction to what we shall calcu- 
late. That would be for p T > 4 GeV/c at y/s = 200 GeV. 
It does not mean that thermal partons are unimportant 
because their recombination with shower partons (the TS 
component) can be dominant (for it production). Both 
TS and SS components are proportional to the probabil- 
ity of hard-parton production, so we may write the basic 
structure of the single-particle distribution in nuclear col- 
lisions in the form 



be denoted by J(k, q, <j>, b) so that we have 

Fi(q, 0, 6) = / dkkfi{k)J{k, q, (f>, 6), (2) 



dN AA 



■(b) = pi(pr,<i>, b) 



pxdprc 



where Fi (q, 0, 6) is the probability that a parton of species 
i with momentum q at azimulthal angle 4> (relative to the 
reaction plane, or more precisely relative to the minor 
axis of the overlap almond in our theoretical consider- 
ation) emerges at the surface of medium. We consider 
only the transverse plane at mid-rapidity. Hi(q,pr) de- 
scribes the hadronization of parton i with momentum q 
to a hadron of whatever type under consideration with 
momentum px- In the recombination model the medium 
effect on hadronization is taken into account by the TS 
component in Hi in addition to the fragmentation term 
that is represented by the SS component [!, 0]- Such 
details are not relevant to our concern here. The main 
point in Eq. (fTJ is that the dependence on <j> and b are 
in Fi(q, <fi, b). The thermal partons can have a minor de- 
pendence of those variables, and is an issue that will be 
discussed later in this paper. 

The essence of our problem is the effect of the medium 
on the hard parton from creation to its emergence on the 
surface. For every hard parton created with momentum 
k the probability of its emergence with momentum q may 



where fi(k) is the distribution for parton momentum k at 
the creation point and has been parametrized in Ref . Q . 
The momentum degradation factor J(k, q, <f>, b) strictly 
should depend on the parton type, but we make the ap- 
proximation of ignoring that dependence, since our focus 
is on the path length dependence. To show explicitly how 
path length enters the problem, we must exhibit the cre- 
ation point at (xq, j/q), which is to be integrated over the 
region of nuclear overlap, weighted by the probability of 
a hard collision. That is, we write 

J(k,q,(f),b) = J dx dyoQ(xo,yo,b)G(k,q,l(x ,yo,<fi,b)) ,(3) 

where Q(xa,yo,b) is the overlap function to be detailed 
below, and G(k,q,£) is the momentum degradation fac- 
tor that changes the parton momentum from k to q in 
a length I. That length depends on the creation point 
(a^Oj^o); the angle 4> that the hard parton is directed, 
and the exit point (xi,yi) for a nuclear collision at im- 
pact parameter b. Again, we defer the calculation of 
£(xo,yo,4>,b) until later so that our description of the 
general scheme is not interrupted. 

Putting Eqs. © and ([3} in (1), we can break up the 
long equation into four parts by introducing J <i£<5(£— 7^?), 
and write 



Fi(q,<f>,b) = 
Fi(q,0 = 



pf+ ss ( PT ,0,6) 



dx dy Q Q(x ,yo,b) 
xS(€-"/£(x ,y ,(t>,b)), (4) 

de^(e,^ 6)^(9,0. ( 5 ) 

dkkfi(k)G(k,q,0, (6) 
/ -£i^,&)^(?,PT).(7) 



A parameter 7 is introduced in Eq. ^ to represent the 
dynamical effect of energy loss on the variable £ that we 
call the dynamical path length. If there is no energy 
loss, then the medium is transparent and £ = what- 
ever the kinematical path length I may be. The value of 
7 will be determined by phenomenology later, after the 
form of momentum degradation, specified by G(k,q,£), 
is discussed. The equations from (J31) to (JTJ) exhibit the 
different parts of the production process separately: (0| 
shows the probability of having a dynamical path length 
£ for a parton directed at (j>; §5§ gives the distribution of 
a parton with momentum q at the surface at angle 0; © 
is the corresponding distribution after traversing an ab- 
sorptive distance £ in the medium; ([TJ is the final hadron 
distribution at px and <f>. A point to note is that the 
relationship between <f> and b that depends on the initial 
configuration is completely contained in Eq. Q. 
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We now define the quantities introduced earlier in Eq. 
(UJ in the framework of the Glauber model for AB col- 
lisions. The thickness function Ta(s) normalized to A 
is 



T A (s) = A J dzp(s, z), J d 2 sT A {s) = A, 



(8) 



where p is the nuclear density normalized to 1 

p{r)=p [l + e^y 5 ]- 1 (9) 

with ro = 6.45 fm and S = 0.55 fm for Au. We shall use 
Ra = 7 fm for the effective nuclear radius. For any point 
(x, y) in the transverse plane, we have 



(x + b/2) 2 + y 2 



(10) 



z\ = 1-s 2 , z l = l-\g-b\ 2 , (11) 



where all lengths are scaled to Ra- The longitudinal 
lengths of A and B at (x, y) are 



2 r z A,B 

L A ,B(x,y) = — dz p(s,z), 

PO J~Z A ,B 



(12) 



where po = 0.285. With a being the inelastic nucleon- 
nucleon cross section, 2a (s) and La(x, y) are related by 

M 

<jT a {s) = uL A {x 7 y), uj = aApo = 4.6 (13) 

for A — 197. The probability for producing a hard par- 
ton at (xq, yo) is proportional to Ta(s + b/2)Tg(s — 6/2), 
where (xo,yo) are the Cartesian coordinates of s, so the 
normalized Q{xo,yo,b) is 

Q(x Q: y ,b) = TA(x 0l y ,~b/2)T B (x 0l y b/2) 

where Tb(xq, yo, 6/2) is the thickness function of nucleus 
B whose center is at an impact parameter 6 on the x axis 
from that of A, located at x — —6/2. 

The path length from the point (xo,J/o) of creation to 
the exit point (xi, y\) on the boundary is 

rti(x ,vo,<f>,b) 

e(x ,y ,<(>,b)= / dtD(x(t),y(t)), (15) 

Jo 

where the integrand is weighted by the local density along 
the trajectory marked by t, since energy loss is propor- 
tional to the medium density, and our aim here is to 
calculate the geometrical part of the factors that con- 
tribute to momentum degradation. Apart from an overall 
normalizaation constant, that density inside the almond- 
shaped overlap region is given by 

D a (x,y) = LoL A {x,y)[l-e^ L ^\ 

+ uL B {x,y)[l-e- wL *<-*»\ (16) 

The normalization is not of concern here because 
(■{xo, yo j <P, 6) appears only in conjunction with 7 in Eq. 



((4]), where 7 is an adjustable parameter to be fixed later 
in fitting the data. 

The exit point (x\,y\) can be determined in terms 
of (xo,yo) an d 4>- The initial point is inside the over- 
lap region, where Q(xo,ya,b) is non- vanishing. The exit 
point is somewhere on the elliptical boundary, since some 
time elapses before the hard parton reaches the surface. 
There is no rigorous way to calculate the transition from 
the almond region to the ellipse, since hydrodynamics 
is not valid before local equilibrium is established; be- 
sides, semihard partons and the ridge formation that they 
give rise to can significantly change the boundary, unac- 
counted for by collective flow [ltj. We determine 
(xi,yi) and the local density D(x,y) along the trajec- 
tory in the following way. First, we take the boundary of 
the ellipse to satisfy 



© 2 + (f) ; 



= 1, 



where 



= 1 - 6/2, h = (1 - 6 2 /4) 



1/2 



(17) 



(18) 



independent of the transit time from (xo,yo) to (x\,y\). 
We map the density D a (x,y) for the almond region to 
the density D(x,y) of the elliptical region by the identi- 
fication 



D(x,y) = D a (xdi/d 2 ,ydi/d 2 ) 



(19) 



along any radial direction measured from the origin 
where x = y = 0. The distances from the origin to 
the almond and elliptical boundaries are denoted by d\ 
and d 2 , respectively, and described in the Appendix. The 
ratio d\jd 2 is very nearly 1 for almost all angles, since 
w and h are the x- and y-axis intercepts of both bound- 
aries. We take the trouble to do the mapping in Eq. (|T9]) 
for the short duration before thermal equilibrium is es- 
tablished in order to render Eq. (|15[) well defined, where 
the upper limit is at the boundary of the ellipse. If the 
point (xo,?/o) is far from the boundary, the transit time 
for the hard parton to arrive at the medium boundary 
can be large so the ellipse would be larger than that de- 
scribed by Eq. (fT7|) . For large transit time the dynamics 
of medium expansion predominantly in the longitudinal 
direction should be taken into consideration. In that case 
we rely on the dynamical scaling behavior found in Ref. 
[llT | to regard our procedure as being insensitive to that 
expansion. Specifically, it means that if the RHS of Eq. 
(jTTJ) were larger, d 2 would be larger, and the rescaling of 
D(x, y) would redistribute the density accordingly and 
the net result of the integration in Eq. (fl~5|) would not 
deviate by too much. D a (x a ,y a ) in Eq. (Til))) is identified 
with the density given in (ITB")) for every point (x a ,y a ) 
inside the almond. By restricting Eq. (| 15[) to only the 
short time duration when (|19p is valid, it is much easier 
to calculate the exit point (xi, y\), where the hard parton 
trajectory along cf> intersects the ellipse, i.e., 



Xi = Xq + t\ cos ( 



Vi = Vo + h sin ( 



(20) 
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where t\(xo, yo, 4>, b) can be determined by solving Eq. 
(|17[) . as done in the Appendix. That is what we use as the 
upper limit of the integration in Eq. (|15j) . thereby giving 
dependence onso, yo, <t> and b to the length 1(xq, yo, 4>, b). 
The result gives an estimate of the path length that 
may possess more general validity than what the static 
approximation of ignoring explicit transverse expansion 
may imply. In this way 4>, b), as defined in Eq. ([4"|l. 
is completely calculable from geometrical considerations. 



III. SCALING BEHAVIOR 

We now focus on P(<£, </>, b) and try to extract as much 
information as possible about its dependence on cj> and b 
before proceeding to the dynamical problem of momen- 
tum degradation and particle production. From Eqs. Q 
and (|14p P(£, 4>, b) is properly normalized as 



J d£P(t,<f> 7 b) = l, 



so the mean £ is 



£(</>,&) = / d^P^^b). 



(21) 



(22) 



Since centrality is more directly accessible in experi- 
ments, we shall freely change from the dependence on b to 
that on c, which is our symbol for centrality; for example, 
c = 0.05 stands for 0-10% centrality. The relationship be- 
tween b and c is well established, and is tabulated in Ref. 
fT3 |. for instance. 
Using Eq. (Q} we obtain 

f(0,c(i>))=7 J dx dyo£(x ,yo,cj),b)Q(xo,yo,b). (23) 

In the next section we shall show that 7 = 0.11. For our 
general discussion here 7 need not be specified, although 
the numerical value will be used when plots are presented 
in figures. We show in Fig. 1 the dependence of £ on <p 
for six values of c for AuAu collisions at ,/s = 200 GeV. 
Except for a shift in the magnitude, that dependence 
seems to be insensitive to c for c ranging from 0.05 to 0.55, 
which corresponds to b from 0.48 to 1.58 in units of It 
suggests that there is a universality in the <p dependence 
for c > 0.05. That universality cannot be exact, since we 
know that when 6 = there can be no dependence on 
<f>. However, for c > 0.05 approximate universality seems 
valid in Fig. 1. If so, then the question becomes whether 
the probability distribution P(£, </>, b) that is a function 
of three variables can more economically be expressed in 
terms of few variables. 

Our first step in addressing that question is to ask 
whether, at fixed b, the distribution is a scaling function. 
More specifically, we define 
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FIG. 1: Average dynamical path length £ vs <j> for six values 
of centrality. 



and ask whether P(£, <fr, b) can be written in terms of a 
scaling function ^{z) in the form 



P(t,4>,b) = ij(z)/Z(<l>,b) 



such that 



dzip{z) 



dzzip{z) 



(25) 



(26) 



(24) 



Distributions that have such properties are referred to 
as satisfyin g K NO scaling, well-known in multiparticle 
production jl3l . ITij . We show in Fig. 2 the behavior of 
ip(z) for (a) c = 0.05 and (b) c = 0.55; in each case there 
are six values of 0: <j> — n7r/24, n = 1, 3, 5, • • • 11. We use 
c for classification, instead of b, because the data to be 
shown in the next section are given in terms of centrality. 
What is notable is that in each case all dependencies on 
<f> collapse to one universal curve in terms of z, when 
expressed in the format of Eq. (|25|) . That is scaling in <f>. 
The next question is whether there is scaling in c. 

To quantify the dependence on c, we fit il>(z) by a 2- 
parameter formula 

V(z) = C ai (l ~ C) a2 /B( ai + 1,02 + 1), C = z/2A, (27) 

where B(a,f3) is the Beta function. The best fits yield 

(a) c = 0.05, oi=Q.37, a 2 = 0.81, (28) 
(6) c = 0.55, ai = 0.57, o 2 = 1.05. (29) 

They are shown by the (red) dashed lines in Fig. 2(a) and 
(b). The (green) dashed line in Fig. 2(b) is a reproduction 
of the (red) dashed line in Fig. 2(a) for the purpose of 
comparison. Visually, the two scaling curves do not differ 
by too much, although numerically the values of a\ and 
ci2 in Eqs. (|28|) and (|29|) are noticeably different between 
the two cases. To a 10% accuracy, one may regard ip{z) 
to be universal for c ranging from 0.05 to 0.55 and <f> from 
to 7r/2. That offers a remarkable degree of simplicity 
in the complex geometrical problem of nuclear collisions. 
In words it means that P(£, 4>, c) depends only on £ and 
£, not on (f> and c separately. 
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(a) c=0.05 
<tH=rm/24, n=1,3,...,11 




FIG. 2: (Color online) Scaling function ip(z) for six values of 
<j> at (a) 0-10% centrality and (b) 50-60% centrality. The red 
lines are fits of the two figures, parametrized by Eq. (|28|) and 
(|29|) ; the green line in (b) is a reproduction of the red line in 
(a). 



IV. NUCLEAR MODIFICATION FACTOR 

Having found some general properties of the collision 
system that are basically geometrical, we now proceed to 
the problem of hadron production and see how the geo- 
metrical insight gained in the preceding section can help 
us organize the dependence on 4> and c. As a consequence 
we shall find an efficient way to compare our result with 
the data on Raa{pt, 4>i c )- 

Let us start by inserting some details that are omitted 
in Sec. 2. Since the PHENIX data on Raa are for 
7T production, we review the formalism for calculating 
the single-pion inclusive distribution 3]. We restrict our 
consideration to midrapidity and write the invariant dis- 
tribution in the recombination model for a pion in the 
direction p in the 1-dimensional form as 



\dN«r 



P 



dqi dq 2 



Fqg( ( 1^92)Rw(qi,q2,p), (30) 



dp J qi qi 
where the qq distribution is in general 

F q q{ qi ,q 2 ) = TT + TS + SS, 

and the recombination function (RF) 



Rn(qi,q2,p) 



(31) 



(32) 



To reproduce the exponential behavior of the observed 
distribution at low pr by TT recombination, we have 
assumed the thermal parton distribution to have the form 



T{qi) = qi 



dN* h 



dqi 



Cqie 



-91 IT 



so that the thermal pion distribution is 
dN^_ 

6 



r 2 

—e~ PT ' T . 



p T dp T 

For the shower parton distribution we have 
S(q 2 ) = J *£*i(ff)^(*>M 



(33) 



(34) 



(35) 



where Sj (z 2 ) is the matrix of shower parton determined 
from the fragmentation functions [7j, and Fi(q) is the 
hard parton distribution at the medium surface before 
fragmentation. The TS contribution to the final pion 
distribution is 



dNj s 
p T dp T 

TS(g,p T ) 



T-E [ -Fi(q)fS(q,p T ), (36) 
Pt j J 1 

^S{(q 2 /q) [ dq lC e-^ T 

qi J 
xR(q 1 ,q 2 ,p T ). (37) 



Finally, the SS component that is equivalent to the frag- 
mentation component is 

A: " ' x ' ' ; " — ' ' (38) 



PrdpT 



Pt, J 1 1 \1 J 



where D{z) is the fragmentation function. When these 
equations are compared to Eq. ([T]), it should be clear 
what Hi(q 7 pT) in the latter stands for. 

The above is a summary of the basic formulas when </> if 
averaged over. Now, as we consider 4> and b dependencies 
explicitly, let us revisit the recombination process. A 
quark q and an anitquark q with momenta q\ and q\ are 
unlikely to recombine to form a pion, if q\ and q 2 are not 
approximately parallel, since their relative momentum 
normal to the pion momentum should not be larger than 
the inverse size of the pion. Thus we shall continue to use 
the 1-dimensional description of the RF. However, the <f> 
distributions of the coalescing partons should be taken 
into account. 

For the thermal partons a description is given in Ref. 
2] that can simultaneously reproduce three empirical 
facts, which are: (a) Raa{pt, <P, c) at 1 < pt < 2 GeV/c 
0, (b) v 2 (p T ,c) [if, and (c) ridge yield Y R {(j) s ,c) 
The physics developed there involves a consideration of 
ridge formation with or without triggers. We omit a de- 
scription of that here, since the details are not of crucial 
importance to our calculation below. We simply adopt 
the thermal parton distribution in the form 

7W, 6) = C(b) qi [l + \aD{x', y', 6)5(0, &)] e ^/ T (39) 
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where a small ^-dependent term is included in the square 
brackets that represents the effect of semihard scattering. 
D(x',y',b) is as defined in Eq. (fT6|) with (x',y') being 
evaluated at 0.17 from the boundary of the almond over- 
lap on the x axis. S((f>, b) is a surface factor defined by 



S(<f>,b) = h[E(0 2 ,a)-E(fi 1 ,a)], 



(40) 



where E(9, a) is the elliptic integral of the second kind, 
9 t = tan _1 (^ tan </>;), </>i = <j> — a, fa = <j> + a, a = 0.33, 
and a(b) = 1 — w 2 /h 2 . It causes a small enhancement of 
the thermal partons due to semihard partons with T = 
0.317 GeV and a = 0.47 in Eq. J3UJ). C(b) is given in 
terms of iV pa rt as 



C(N paxt ) = 1.1N%* GeV" 



(41) 



These complications are included here for completeness, 
even though their effect on the final result is small. 

The shower partons, on the other hand, depend cru- 
cially on </>, as a comparison between Eqs. ([7]) and (|35|) 
reveals the presence of 4> in Fi(q, <f>, b) in the former, which 
is the crux of the problem at hand. In Eq. ([2|) the dis- 
tribution of a parton with momentum q at the medium 
surface is related to the parton distribution fi(k) at the 
point of creation whose dependence on k is known Q. 
The relationship between fi{k) and Fi(q,<p,b) is compli- 
cated in view of Eqs. © and ([3]), but has been simplified 
to Eq. ([6]) by the introduction of £. In terms of £, the 
degradation of momentum from k to q can be written as 
a simple exponential for 6 < k < 12 GeV/c, i.e., 



G(k,q,£) = q5(q-ke-t), 



(42) 



whose justification based on an approximation of QCD 
result is given in Ref . @ . Using this in Eq. © yields the 
simple connection 



E(q,0=q 2 ^Mq^) 



(43) 



With this substituted into Eqs. (j3"6")l and (pi5)) in place of 
Ei(q), we obtain the TS and SS contributions to the pion 
distribution as functions of px and £ 



Pi 



TS+SS 



(Pi 



,0 = JjJ2 F ^® H ^>Pt)' ( 44 ) 



where Hi(q,pT) is as noted after Eq. (|3"5)l . To deter- 
mine the dependence on 4> and b, it is necessary to make 
one last transform according to Eq. ([5]), using what we 
have learned about P(£,(j),b). It should now be clear 
that whereas Eqs. (|4"3"]l and (|44l) describe the dynamical 
processes of hard parton scattering, momentum degra- 
dation and hadronization, the geometrical complication 
is contained in P(£, 4>, b). The validity of this separation 
remains to be verified by comparison with data in a way 
that can best reveal the role of £. 

We calculate the nuclear modification factor 



Raa(pt,4>,c) 



dNl A /d PT dc^ 
N con dN£ p /dp T 



(45) 



where N co u is the average number of binary collisions and 
dNp p /dpx is the pion distribution in pp collisions. For 
AuAu collisions we include the TT component, although 
it is small compared to the TS and SS components for 
Pt > 4 GeV/c. In Ref. 1] good quality data are given 
for 4 < pt < 5 GeV/c, six bins of 4> and five bins of 
c. At higher px the errors are larger. We calculate the 
same quantities for px = 4.5 and 7.5 GeV/c, and for 
similar ranges of <fi and c. An important step we take is 
to present our results as functions of £(</>, c). 
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FIG. 3: Theoretical result showing scaling behavior of Raa 
at pt = 4.5 GeV/c for six values of centrality. The points for 
different <f> values between and 7r/2 have the same symbol 
when the centralities are the same, with <f> increasing from left 
to right with n = 1,3, •••,11. The solid line is a fit by Eq. 
(|47p and the dashed line is an exponential fit by exp(— 2.6f). 

Since for each pair of values of <j) and c, we can cal- 
culate R aa {pt, 4>, c) and £(0, c), they can therefore be 
plotted one versus the other, as in Fig. 3. We have done 
this for six values of (f> at n7r/24, n — 1, 3, • • • , 11, and 
six values of c at 0.05m, m = 1, 3, • • • , 11, with pt fixed 
at 4.5 GeV/c. Points for various 4> with the same c have 
the same symbol; they are spread out from left to right 
with increasing values of <p. All 36 points lie remarkably 
well on one universal curve. It means that different ex- 
perimental conditions on <j> and c yield the same value 
of R\a so long as their values of £,(<$>, c) are the same. 
The nature of the behavior of R\ A (0 depends on the 
details of the dynamics. It can be approximated by an 
exponential 



Raa(0 



Pt=4.5 



exp(-2.60 , 



(46) 



as shown by the dashed straight line in Fig. 3, although 
a better fit is 



Kaa(0 |pt=4.5 



1.22(1 + ^/0.69) 



-2.86 



(47) 



shown by the solid line. That is a simple formula for a 
wide range of <fr and c. 

Similar calculation can be done for pr = 7.5 GeV/c. 
Fragmentation dominates at such a high px- We include 
all contributions to the hadronization with the result 
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FIG. 4: Same as for Fig. 3 but for p T = 7.5 GeV/c. The solid 
line is a fit by Eq. 
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FIG. 5: Raa vs iVpart for 3 values of <f) at pT~i-5 GeV/c. 
The data are from Ref. Q. The lines are from calculation in 
the recombination model. 



shown in Fig. 4. Again, good scaling is found for all 
(j> and c. The solid line shows a fit of the £ dependence 
by 



Raa(Q |pt=7.5 



1.25(1 + ^/0.58) 



-2.62 



(48) 



A comparison with Eq. (|47|) indicates that there is no sig- 
nificant change in R\a(0 as pt is increased from 4.5 to 
7.5 GeV/c. Average values of the two sets of parameters 
can yield acceptable fits of all points calculated. Thus we 
learn that there is an approximately exponential suppres- 
sion of Raa in terms of the path length £ for all <\> and 
all centrality < 60%, for pt > 4 GeV/c, bearing in mind 
that £ is no greater than 0.65. The degree of suppression 
is no more than a factor of 5, which is a familiar fact. 
In Figs. 3 and 4, we see clearly how that suppression is 
related to one another at different ci> and c. 

To compare our calculations with the data we have one 
basic parameter to vary, which is 7 in Eq. Q that relates 
f to £(xo, 2/o,0, b), defined in Eq. ([IS]). That integral for 
£(xo,yo,4>,b) is along the path of the hard-parton tra- 
jectory weighted by the local density D(x,y), which is 
proportional to the longitudinal length at (x,y). With 
the normalization factor being absorbed by 7, £ itself 
has the characteristic of length, and 7 behaves as inverse 
length, though all length are in units of Ra- The dimen- 
sionless £ = j£ is the dynamical path length, in terms 
of which the momentum degradation can be expressed 
as e~£. Our formalism does not provide a way to cal- 
culate 7, but it does account for the geometrical details. 
In QCD energy loss can be calculated, but without the 
details of the medium. If the scaling behavior that we 
have found can also be found in the data, then the value 
of 7 can be extracted, giving justification to our proce- 
dure of condensing all the geometrical complications into 

In Fig. 5 we show the data from Ref. |l| on Raa vs 
A part at 4 < pt < 5 GeV/c for 3 values of </>. We adjust 
the parameter 7 to fit the data and obtain 



There is another parameter that has not been mentioned 
so far. It is the lower limit to of the integral in Eq. ([TJJ. It 
can differ slightly form to take into account the absence 
of immediate effect of energy loss at the very beginning of 
a hard-parton trajectory. A non-zero to shortens the ef- 
fective path length and can increase Raa at small -ZV part 
and small cf>. A better fit of the data in that region is 
achieved by using to — 0.025, which corresponds to less 
than 0.2 fm, a value small enough so that it does not play 
a significant role in the general description of the prob- 
lem. All figures shown above are results of calculations 
done with 7 and to given these values. Another way of 
presenting the results is to plot Raa versus 4> f° r differ- 
ent centralities. That is done in Fig. 6, showing good 
agreement with data 
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FIG. 6: Raa vs <f> for 5 bins of centrality at pt=4-5 GeV/c. 
The data are from Ref. Q. The lines are from calculation in 
the recombination model. 

As we can see from these two figures, the data plotted 
in terms of </> and c vary over wide ranges of values, ex- 
hibiting no organized simplicity. It is therefore important 
to replot Raa versus £(4>, c) to check the scaling behav- 
ior found in our theoretical study. Fig. 7 shows the data 
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replotted in that format; the solid line is a direct trans- 
fer from the solid line in Fig. 3 that fits the theoretical 
points by Eq. (|47|) . It is a way to show the relationship 
between the two figures, but in actuality the theoretical 
and experimental points agree with one another better 
than how the analytic formula can represent them. Given 
the experimental errors, the scaling behavior is evidently 
in the data. 
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FIG. 7: Experimental data on Raa for pr=4-5 GeV/c from 
Ref. [l[ are plotted against the scaling variable £. The (j> 
values are the same as described in Fig. 3. The solid line is 
the same as the solid line in Fig. 3 that represents the result 
from theoretical calculation. 
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FIG. 8: Experimental data on Raa for pr=7-8 GeV/c from 
Ref. [l| are plotted against the scaling variable £. The (j> values 
are the same as described in Fig. 3. The solid line is the same 
as the solid line in Fig. 4 from theoretical calculation. 

For pt in the 7-8 GeV/c range we show the data only 
in the scaling format, as in Fig. 8. The solid line is a 
reproduction of the solid line in Fig. 4, described by Eq. 
(|48|) . Within errors which are large, the agreement be- 
tween theory and experiment is good in both the magni- 
tude and the £ dependence. No parameters have been 
adjusted. Thus we learn not only that scaling is ro- 
bust, but also that the dependence on £ changes with 
Pt mainly because the SS component of hadronization 
becomes most dominant at higher px , since no other fac- 
tor besides Hi{q,pT) in Eq. (|44")) depends on pt- The 



medium effect is primarily contained in Fi(q,£). 

It should be mentioned that the path length depen- 
dence has also been investigated in Ref. [1]. Scaling be- 
havior has been found in different variables for different 
Pt ranges. Their definition of pL xy / p ccnt in fm is closest 
in spirit to our dimensionless £. The linear dependence 

— 1/2 

of In Raa on the path length L and Uj, has also been 
checked experimentally in Refs. (l7l. fl8j. as predicted in 
QCD @. Since our £(0, c) is far more detailed than L 
and our px range includes the contribution from thermal- 
shower recombination, our result cannot be related sim- 
ply to any QCD formula. 

It is interesting to show how the £ dependence can 
be obtained in a more direct and transparent way, when 
an approximation is made at high pt- The hard-parton 
distribution fi(k) has been parametrized in Ref. Q in 
the form 



f t (k)^(l + k/k t )-^ 



(50) 



where for i = g (u),ki = 1.77 (1.46) GeV/c and /?i = 
8.61 (7.68). For simplicity, we approximate the formula 
by a simple power law for k > 7 GeV/c 



fi(k) oc k-K, 



(51) 



where f3' g = 7.8 and /3' u = 7.08. The advantage is that 
becomes, from Eq. (|43]) , 



= # - 2, 



(52) 



which renders the q and £ dependencies separable. We 
can then perform the q integration in Eq. (I44[) indepen- 
dent of £, and obtain the gluon and u-quark contributions 
separately as 



RUa(Pt,(I>,c)<x J dSP&hje-f* 



(53) 



The scaling property of P(^,cj),c), as described by Eq. 
(|25| . can now be used to yield 



r aa(pt,0 oc / dztp(z)e 



(54) 



which can simply be evaluated in conjunction with Eq. 
(|2T|) . The result is shown in Fig. 9 for gluon in (green) 
solid line and for u quark in (red) dashed line, the nor- 
malizations for which are freely adjusted. It should be 
noted that the curve in Fig. 7 is from complete calculation 
taking the contributions from all partons into account 
without using explicitly the scaling property of -P(£, (j), c), 
while the two curves in Fig. 9 are obtained through the 
approximation of Eq. (|51[) so as to exhibit directly the 
roles of (a) the parton distribution function fi(k), (b) 
the exponential relationship between k and q, and (c) 
the scaling function ip(z). All three of these properties 
are the main content of the integrand in Eq. (|54l) that 
shows clearly the three basic aspects of the physics that 
give rise to the observed behavior of Raa- 
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FIG. 9: (Color online) The same data as shown in Fig. 8 for 
Pt=7-8 GeV/c. The (green) solid line is for gluon fragmenta- 
tion only and the (red) dashed line is for M-quark fragmenta- 
tion; the normalizations for both contributions are arbitrarily 
adjusted so as to show the £ dependence compared to data. 



It is of interest to point out the implication of the two 
factors in the integrand of Eq. (f53")l . or (|54l) . The factor 
exp(— £ri) is basically an expression of the hard-parton 
distribution fi(k) in terms of £, indicating that the prob- 
ability of a hard parton contributing at large £ to an ob- 
served pion is exponentially suppressed. <j>, c) gives 
the probability of creating a hard parton at £ from the 
surface. Their product, expressed as ip(z) exp(— z^r.;) in 
Eq. (l54l) . shows that, although the mean z is 1, the small 
z portion dominates due to higher suppression at large 
z. It means that most of the partons that contribute 
to pion production are created near the surface. If the 
result of the integration over z can be approximated by 
exp(— z^Ti) with z = 0.5, we see that the dependence of 
Raa(0 011 £ i s roughly the exponential behavior given 
by Eq. (gBJ). 



V. BACK-TO-BACK DUETS 

Having seen how the dynamical path length £ plays a 
crucial role in determining the general behavior of Raa, 
it is natural to extend the study to consider back-to- 
back dijets, whose properties clearly depend on the dis- 
tances that the two generating hard partons traverse in 
the medium. That problem has been studied previously 
in Ref. [f| without taking into account azimuthal depen- 
dence explicitly. We now can investigate the effect of 
£ scaling and find a more transparent way to show the 
dependence of away-side yield per trigger on <fi and c. 

Using the same notation as in [4] where pt and pb de- 
note the momenta of the trigger and the associated par- 
ticle on the away side, respectively, at mid-rapidity, we 
have by simple generalization from Eq. (JTJ) 



dN AA 



ptpbdptdpbc 



-(b) = P2(pt,Pb,<f>,b) 



/" — — ©, 0, b)Hi( qi , q 2 , 4>, 6), (55) 



where q\ is the parton momentum on the near-side sur- 
face, and q 2 on the away-side surface. They are re- 
lated to the momentum k at the creation point by 
G(k,qj,£j), j = 1,2. Eq. (j4]) is now generalized to 



P{tii,&,<l>,b) = J dx dy Q(x ,y ,b) 

x5(£i - 74 (x , yo,<f>, b))S(& - 7^(^0,2/0,0, ^X 56 ) 

where £i(xq, yo, </>, b) is as defined in Eq. (|15l) . and 
•^2(^0, 2/o, </>, b) is defined similarly but in the opposite di- 
rection with the upper limit of integration replaced by 



h{x ,y ,4>,b) =t 1 (x ,yo i <j> + 'K,b). 



(57) 



In our calculation below we shall set the lower limit of 
both integrals at t = 0, ignoring the small to mentioned 
after Eq. (|I9"|). As in Eqs. © and @ we have 

F i {q 1 ,q 2 ,cj>,b) = J dfrdfri^i, 6, ^b)F l (q 1 ,q 2 , ft, £ 2 )(58) 

Fi(qi, q 2 , Ci, 6)= / dkkf i (k)G(k,q 1 ,Zi)G(k,q 2 ,bl59) 

Computation of pi(jPtiPbi4>ib) can now be performed 
from the equations given above without questioning the 
scaling properties. However, we have learned from the 
preceding section that considerable clarity can be gained 
by presenting the result in terms of the scaling variable 
I 

Our next step in that direction is therefore to study 
the properties of P(£i, £2, <f>, b). We define 



^•=0/£(<M), J = 1,2, 



(60) 



where £(0, b) is the same as determined in Sec. 3 and 
shown in Fig. 1. It should be noted that although 
Eq. (|2"3"]l expresses £(</>,&) in terms of the near-side path 
length £i(xo,yo, 4>i b), it is the same as for the away side, 
since (xo,2/o) is integrated over the entire overlap so re- 
placing <f> by 4> + T does not lead to any difference. We 
now ask whether -P(£i, £2, 0, b) has the scaling behavior 



where 



P{liiA2A,b) = ^{zx,z 2 )/?{<t>,bl 



dz 2 ^(z 1 ,z 2 ) = ipizip). 



(61) 



(62) 



For (j> — 7r/24 and c = 0.05, we have £ = 0.55 and we 
can plot -P(£i,£2) in the format of ^>(zi,z 2 ) vs (zi,z 2 ), 
as shown in Fig. 10. It is essentially symmetric in z\ 
and Z2, as it should, and peaks near the boundary of a 
maximum z\ + z 2 . The existence of such a boundary is 
due to the finite size of the medium that puts an upper 
limit on £1 + £2. Fig. 11 shows the maximum z 2 as a 
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function of z\. Evidently, the symmetry in z\ and z<i 
suggests the use of the variables 

z± = z 1 ±z 2 , (63) 



.max 



= 2.4. We show in Fig. 12 the projections 
of ^(z+jZ-) on z + for 4 values of Z-. The four curves 
are in essence very similar although in detail they cannot 
be identical because z\ and z-i must both be positive so 
z + has different lower limits depending on z_. 
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FIG. 10: (Color online) A 3D plot of the scaling function 
ty(zi,z 2 ) for c = 0.05 and (j> = 7r/24. 
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FIG. 11: Maximum z 2 vs z\ for c = 0.05 and 4> = 7r/24. 

Figures 10-12 are for <p = tt/24 and c = 0.05. To 
show scaling it is necessary to do the calculation for other 
values of cf> and c and demonstrate universality. We have 
done that and found plots like Fig. 12 that have almost 
no perceptible differences among them and are thus not 
exhibited. We conclude that P(£i, £2, <t>, b) does satisfy 
scaling as described by Eq. (IBTj) . 

We can now return to the dijet spectra and present the 
results on yield per trigger defined as 



PikPuVbA-, b) 
pi(p t ,(j>,b) 



(64) 



which is plotted as a function of £ in Fig. 13 for p t =8 
GeV/c and pb = 4,6,8 GeV/c. The symbols are the 
same as those in Fig. 3, representing (f> = n7r/24 and 
c = 0.05m. There is clearly good scaling behavior in £. 



FIG. 12: The scaling function ^(z+, z~) vs z+ for four values 
of z- at c = 0.05 and <j> = n/24. 



They arc the results of direct calculation of P2 ijpt iPbi^i c ) 
without using the properties of ^(zi, 22) m Figs. 10 and 
12. The dependence on £ is similar among the three sets 
of pi, values, although for pb — 4 GeV/c the decrease with 
£ is slightly slower than for 6 and 8 GeV/c. The normal- 
ization of each set is due to many factors, among which 
are obviously the fragmentation functions Di(pb/q2) on 
the away side contained in Hi{qi 1 q2,pt,Pb) in Eq. (|55l) ; 
they can directly account for the decrease with increas- 
ing pi, if q2 were held fixed. Since £2 affects the depen- 
dence on q2, it is understandable why the £ dependence 
of Y^ ay (p t ,pb,0 at p b — 4 GeV/c can differ from that 
at pb — 8 GeV/c, though only slightly. 
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FIG. 13: Away-side yield per trigger at trigger momentum 
pt = 8 GeV/c and for three values of the associated- particle 
momentum p^. The symbols are the same as in Fig. 3 for 
various values of c and if). 

Another way of presenting the results shown in Fig. 13 
is to exhibit the pb dependence for fixed values of £. In 
order to compare with data we relate Y a ™ ay (p t ,pb) to the 
trigger-normalized fragmentation function D^^zr) p^j . 
where z T =p b /p t , by 



a** iz T , = pmY^T ( Pt , Pb , , 



(65) 
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so that the integral of both sides over zt gives the same 
total yield of pions on the away-side of a fixed trigger 
pion at pt,<fr,c. In Fig. 14 we show the result for three 
values of £. The Zt dependence does not change very 
much in the range 0.5 < zt < 1.0 corresponding to 
4 < pb < pt = 8 GcV/c, when £ is increased from 0.1 
to 0.5. There are recent data from STAR on 7r°- cor- 
relation on the away side for 8 < pt < 16 GeV/c and 
0.35 < z T < 0.95 at 0-10% centrality 0. They are 
shown in Fig. 14 for qualitative comparison with our cal- 
culated result, although the large variation in trigger mo- 
mentum may render quantitative comparison meaning- 
less. The general zt behaviors in theory and experiment 
are consistent. 
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Writing ^(2+, Z- = 0) as $(z+), we then have 
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FIG. 14: Trigger-normalized fragmentation function on the 
away side. Data are from Ref. [2(| 



The value of £(</>, c) is determined from geometrical 
considerations, apart from the multiplicative factor 7. 
The dependence of the yield on £ is, however, rooted in 
the dynamics of hard scattering and energy loss. To shed 
some light on the role played by scaling, we go first to 
Eq. (|59"1) and with the help of Eq. (|4"2")) write it in the form 

i ? i(?i,32,a,6) = [ k3 M k )] k=iqiq2 e^i2 V /2 

x%ie 6 - g 2 e« 2 ). (66) 

Using Eqs. and (joTj) we obtain 



Fi(qi,q2,Q = J dz 1 dz 2 ^(zi,z 2 ) 

x [fc 3 /i(fe)], =(?i?2e | ( , 1+ , 2))1/2 5(gie^- g2 e^).(67) 

For symmetric dijets, pt = Pb = 8 GeV/c, qi and 92 
are on average equal, so z\ ~ z%. Eq. (|6"T)l can then be 
approximated by 

Fi(qi,q2,Z) oc / dz+dz^{z + , z^-^'+^t^Z- 



where Eq. (1511) has been used, and multiplicative fac- 
tors involving q\ and qi are omitted from the expression. 



(69) 



where x I'(z+) is very nearly the solid line in Fig. 12, ren- 
dering the £ dependence explicit. Since *f/(z+) is very 
different from ip{ z )i P2(Pt,Pb, £) is significantly different 
from p\{pt, £), which is of the form given in Eq. (|54p . In 
both cases there is scaling in £ by virtue of the scaling 
properties of P(£, ^, 6) and P(£i, £2, 0, 



VI. TWO-JET RECOMBINATION AT LHC 

In heavy-ion collisions at CERN-LHC the density of 
hard partons created is so high that the shower partons 
from neighboring jets can recombine to form hadrons. 
It is found in Ref. [lj| that large p/n ratio, as high as 
20, can be a signature of such hadronization processes. 
Here we examine a different signal that may show up 
even more dramatically in the preliminary data collected 
in the initial phase of the operation of LHC. In simple 
terms the £ scaling of Raa may be badly broken, and 
Raa itself can be very large. 

We consider pion production only at high px where 
thermal partons can be ignored, but not so high that 
the hard-parton density becomes low. For definiteness, 
we set p T = 10 GeV/c for y/1 = 5.5 TeV in PbPb col- 
lisions. As before, we consider only hadron production 
at midrapidity. The 2-jet contribution to the inclusive 
distribution is a simple generalization of Eq. (jTJ) 



dN 2 A 
prdprdf, 



fdqdq>_^ b)Hu>(q,q',PTm 



where q and q' are the momenta of two hard par- 
tons at the surface of the medium, both directed at 0, 
and Hui(q,q' ,pt) describes how those two hard partons 
hadronize through their separate shower partons that re- 
combine, i.e., 



1 



dqi dq 2 



P T J 9i 92 V 1 J \ 1 j 
xRl( qi ,q 2 , PT ). (71) 

The parton types i and i' are summed over g, it, d, u and 
d, so there are 25 combinations of hard partons, each frag- 
menting into j and j' shower partons. The recombination 
function R^(qi,q 2l pT) must now take into account the 
probability that the two shower partons of types j and j' 
from two jets can coalesce to form a pion. Since the two 
hard partons are parallel, the primary factor is how far 
apart they are. Secondarily, the shower partons in a jet 
may form a cone, so the RF is proportional to the overlap 
of two such neighboring cones. There are complications 



'72 
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that we do not know well enough to describe quantita- 
tively. Fortunately, the main effect that we are searching 
for is independent of the details. We approximate the 
RF by 



(72) 



where i? 7r (q , i, q2,Pr) is given in Eq. (|32l) . and Y is a nu- 
merical factor that we shall allow to have widely different 
values. More specifically, we let 
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If a crude estimate of T is made by considering the dis- 
tance between two parallel trajectories, then it would be 
the ratio of the pion radius to the effective nuclear ra- 
dius, Rn/Rf. Taking R^ to be 0.6 fm and the Rf « 6 
fm (slightly less than Ra due to the predominance of the 
hard scattering point near the center where Q(xo,yo,b) 
is higher), we get V ss 0.1. 

Transforming the (<j>,b) variables into £ as in Eq. ([5]), 
we have 

P?(pr,<t>,b) = J d^P(t4,,b)P(£,',4,,b)p\\p T ^,ei7±) 

where 

J 1 1 77! 



The two variables £ and £' are independent, unlike £i and 
£2 in P(£i, £2, 4>i b) considered in the preceding section. 
But for any <f> and b there is only one £ (0, b) . Thus KNO 
scaling implies 
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FIG. 15: Scaling behavior of Raa at LHC for p T =10 GeV/c 
and 2-jet recombination overlap factor V = 0.001. The sym- 
bols are the same as in Fig. 3 for various values of c and <f). 
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Pi 3 (pr,<j>, b) = / dzdz'il){z)i){z')p 1 J {j>T,z,z',C)^ ( 76 ) 



In our calculation we include also the 1-jet contribution 
in Eq. (|44l) . in addition to the 2-jet contribution above. A 
number of parameters are changed in going from RHIC to 
LHC, such as a™ = 60 mb, A = 207, u = 14, and those 

inci ' ' ' 

in fi (k) [8j , but we have fixed 7 at the value given in Eq. 
(14^1) as a reasonable working hypothesis. To calculate 
Kaa(pt,£) we use dN^/p T dp T = 8.7 x 10- 7 (GeV/c)- 2 
for pt — 10 GeV/c at y/s = 5.5 TeV, getting the result 
shown in Fig. 15 for T = 10 -3 . The symbols are the 
same as those in Figs. 3 and 4. It is clear that there is 
excellent scaling as before. However, the reason for the 
scaling behavior is the dominance of 1-jet contribution. 

To show the effect of the 2-jet contribution we repeated 
the calculation for T — 0.1. The result is shown in Fig. 16. 
Evidently, scaling is badly broken. The root cause can 
be traced to the fact that two jets involve two separate 
creation points of hard partons, the probability for each 
being proportional to N co i\. Thus the N% a dependence of 
Pi (Pt> £> C) is n °t scaled out by the definition of R AA in 
Eq. (|45p . The admixture of and p 2 ^ therefore depends 
implicitly on N co \\ that ruins the £ scaling. To confirm 



FIG. 16: Non-scaling behavior of Raa at LHC for pt = 10 
GeV/c and F = 0.1. The symbols are the same as in Fig. 3 
for various values of c and d>. 



our rationale, we have left out p^ and studied P\ /N^ Q ^ 
from Eq. ([To]) ; indeed, we found scaling. 

The contrast between Figs. 15 and 16 is striking, the 
only difference in the input being the values of Y used. 
Experimental data should readily be able to distinguish 
the two scenarios, and give hints on the value of Y that 
we know too little to calculate. At pr much higher than 
10 GeV/c parton momenta k and k' are higher, so the 
parton density would be lower and the probability of 2- 
jet recombination would be suppressed. At lower px the 
ratio of 2j-to-lj contributions in Pb-Pb collisions without 
the r factor increases by a factor of around 30% when 
Pt is lowered from 10 to 6 GeV/c. Hence, if Y = 10~ 3 , £ 
scaling would persist, but if Y = 10 _1 , then the scaling 
violation would also continue to prevail. 

Another aspect about Fig. 16 that is perhaps more 
spectacular than the scaling violation is the magnitude 
of Raa that can exceed 1 for small cf> (left side of each 
type of symbols). For central collisions (c = 0.05, open 
stars) for which £ is between 0.5 and 0.7, the magnitude 
of R\a w ith dominance by 2-jet recombination for Y = 
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0.1 can be four times larger than the 1-jet contribution 
for T = 10~ 3 , both being for pion at px = 10 GeV/c. 
The increase is not so large for mid-central collisions at 
c = 0.55. Clearly, the N^ oll dependence of boosts R\ A 
in central collisions to completely overturn the effect of 
nuclear suppression. Thus the detection of large R\ A 
(and even larger R AA according to Ref. [Hj]) at LHC is 
an immediate signal of two-jet recombination. 



VII. CONCLUSION 

For hadron px large enough so that the dependence 
of R AA on centrality c and azimuthal angle (f> is closely 
related to the propagation of hard partons in a dense 
medium, our treatment of the geometrical aspect of the 
problem has led to the finding of a scaling property that 
not only simplifies considerably the presentation of such 
dependence, but also offers a clearer picture from which 
one can learn about the nature of medium effect on jets. 
In terms of the dynamical path length £ (0, c) there is 
in R AA an approximate exponential behavior for all <f> 
and c. Based on the results of our calculation in the re- 
combination model, we present the data of Ref. [l| in a 
similar format and found the scaling behavior for px at 
4-5 and 7-8 GeV/c over wide ranges of <p and c. Further- 
more, the dependence on £(0, c) is as we have obtained 
theoretically. 

Experimental data presented in the form shown in 
Figs. 7 and 8 can be done for any pt without theoret- 
ical input. The normalization for £(0, c), as defined in 
Eq. ([23]) . does involve a fitted parameter 7 that is not 
directly accessible from the data. Leaving it out, one 
can calculate the mean path length £((f>, c) from geome- 
try and the Glauber model; thus experimental plots of 
R AA vs I can be made without changing the character of 
scaling. The nature of the I behavior in such a plot does 
depend on dynamical details at the pt considered, such 
as the dominance of TS or SS component in hadroniza- 
tion. Our general conclusion is that for pion production 
R AA behaves approximately as e -2 ' 6 ^ with 7 0.1. This 
is a very succinct statement of what is observed at RHIC, 
but can be grossly violated at LHC. 

Hadron correlation in back-to-back dijets is sensitive 
to the path length of the associated particle on the away 
side. The double-hump structure observed on the two 
sides of Acj) = it [2ll . l22j has recently been examined 
as a function of the trigger 4>, revealing differences in 
the "head" and "shoulder" regions, as well as asymme- 
try in A(f) that depends on path length [23| . The inter- 
esting phenomena unfortunately cannot be compared to 
the findings in our study here because the trigger mo- 
mentum is only in the 2-3 GeV/c range, which is much 
lower than what we consider. At lower px the effect of 
semihard scattering is important, since such scattering 
is pervasive and can influence the medium in ways not 
accounted for by hydrodynamics. Our investigation has 
been on the opposite effect — that of the medium on 



high pr jets. For pt = 8 GeV/c we find that the yield 
per trigger on the away side has essentially a universal £ 
dependence for pi, > 4 GeV/ c. Only the magnitude of the 
yield depends on p\> because of the fragmentation func- 
tions, but the dependence on cj> and c is well described by 
the path-length consideration through £ scaling. When 
expressed in terms of D{zx), the zx dependence seems to 
be in accord with the data. Verification of the £ scaling 
in the head region at high p t and pb would give even more 
empirical support to our findings in dijet correlation. 

At LHC that £ scaling is broken when 2-jet recombi- 
nation becomes important at high hard-parton density. 
If some aspect of our prediction is to be found valid in 
the data, the striking phenomenon to be seen in the pre- 
liminary result from LHC will be large R AA at px ~ 10 
GeV/c. It would be a clear experimental indication that 
the conventional understanding about the effect of dense 
medium on hard partons (i.e., R AA < 1) is inadequate, 
and that the hadronization process through recombina- 
tion of jet showers can obliterate the effect of momentum 
degradation, resulting in R AA > 1. The non-obervance 
of a large value of R AA in excess of 1 does not falsify 
our consideration, since it would be due to the smallness 
of the overlap probability T which we have no reliable 
way to calculate. It would reduce the impact of 2-jet 
recombination, but contributes to better understanding 
of the hadronization of partons at high density. How- 
ever, if the large R AA is seen in the data at px ~ 10 
GeV/c, it would mean that the background of the much 
higher-j»T jets is significantly larger than what can be 
expected from hydrodynamical predictions by extending 
similar considerations from the RHIC regime. 
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Appendix A: GEOMETRICAL DETAILS 

In mapping from the initial density in the almond- 
shaped overlap region to the density in the initial ellipse, 
we have Eq. (TK)1) , in which d\ and d,2 refer to the distances 
from the origin to the almond and elliptical boundaries, 
respectively. They are defined, for any radial <f> r from 
the origin (not to be confused with the parton's <f> from 
xo,Vo), by 

(^ + d% cos0 r ) 2 + (efi sin</> r ) 2 = 1, (Al) 

(— cos0 r ) 2 + (^sin0 r ) 2 = 1. (A2) 
w h 
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The above are for — 7r/2 < cj) r < ir/2; the case for the left 
half plane can be obtained by symmetry. Solving them, 
we obtain 



di 
d 2 



1/2 



"I 1/2 



h 2 



(A3) 



Note that di/d 2 = 1 when <fi r — or tt/2, as they should, 
since w and h are the x and y intercepts of the almond as 
well as the ellipse. For general <p r , d\j di is not far from 1, 
but the small difference makes the density D(x(t),y(t)) 
in Eq. (|15[) well defined for any (xq, j/o) inside the almond 
and the corresponding ii(xo, 2/o, 4>i b) on the ellipse. 

To determine ti(xo,yo,4>,b), we substitute Eq. ([20)1 
into (|T7|) and solve for t\, getting [To| 



(B 2 



AC) 1/2 -B 



/A, 



(A4) 



where 



A : 







( — sin 




( — cos 6 


)'♦ 




)' 


\ w 









is = — COS ( 

vr 



yo . 

h? sm< 



(7=1- (x /u;) 2 - (yo/M 2 



(A5) 



(A6) 



(A7) 



It then follows from Eq. (fl5|) that £(xo,yo,(f>,b) can be 
calculated entirely from geometry and therefore £((/>, b) 
as well, using Eq. (f23f without the multiplicative factor 



7- 
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